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ABSTRACT 



6 

, ^, ' We have carried out a survey of X-ray emission from stars with giant planets, 

combining both archival and targeted surveys. Over 230 stars have been currently 

>- ■ identified as possessing planets, and roughly a third of these have been detected 

OO 
O 
cn ■ main sequence star systems with detected planets, comparing subsamples of stars 

that have close-in planets with stars that have more distant planets. This analysis 

Q ■ reveals strong evidence that stars with close-in giant planets are on average more 

OO ■ 

O 



in X-rays. We carry out detailed statistical analysis on a volume limited sample of 



X-ray active by a factor k, 4 than those with planets that are more distant. This 
result persists for various sample selections. We find that even after accounting 
for observational sample bias, a significant residual difference still remains. This 
observational result is consistent with the hypothesis that giant planets in close 
proximity to the primary stars influences the stellar magnetic activity. 

Subject headings: stellar activity - extrasolar planets - magnetospheres - surveys 
- X-rays: general - methods: statistical 



1. Introduction 

After centuries of ignorance on planetary systems beyond our own, we are now in an 
era when stars are being routinely identified as possessing planets. Since the discovery by 
Mayor & Queloz (1995) of a Jupiter-type planet (mass 0.44 Mj) orbiting close in to the 
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G2 star 51 Peg (0.05 AU), almost three hundred new Extrasolar Giant Planets (EGPs) have 
been found|l| 

One of the surprising results so far have been the detection of planets with orbital 
radii, op, much smaller than seen in the Solar System, with values as low as ap ~0.03 AU. 
At such small separations, it is likely that these giant planets will have measurable tidal 
(oc ap^^") or magnetic (oc ap'"^ for large ap) effects on the primary stars (Cuntz, Saar, & 
Musielak 2000; Saar, Cuntz, Shkolnik 2004). Both these effects increase non-linearly for small 
values of ap, leading to potentially dramatic disruptions of the stellar environment for stars 
with close-in giant planets. How these disruptions affect coronal activity remains an open 
question. It is possible that activity will be enhanced because of the magnetic interactions 
between star and planet. Tidal bulges will also have an effect on the stability and energetics 
of the chromosphere. Hitherto, analyses of variables that control stellar activity levels for 
single stars have focused on rotation rate, mass, age, evolutionary state, and to some extent 
metallicity, as factors determining coronal activity. Now we are presented with another 
possibility, viz., the existence of close-in EGPs may also be a controlling parameter. It is 
well known that binary stars with close stellar companions are generally more active than 
single stars (see e.g., Pye et al. 1994). Stars with close-in EGPs may be of the same type, 
though it is unclear whether similar mechanisms of activity enhancements hold at both 
extremes of this family. 

A preliminary search for planet-induced stellar activity enhancement was carried out 
by Bastian, Dulk & Leblanc (2000) in the radio and by Saar & Cuntz (2001) in the optical. 
But due to insufficient sensitivity, neither study succeeded in uncovering clear evidence for 
this phenomenon. However, a sensitive search for activity (Shkolnik et al. 2003; see also 
Cuntz & Shkolnik 2002) detected enhanced stellar emission in the Ca II HK lines from the 
chromosphere of HD 179949, in phase with the orbit of its close-in giant planet (Porb = 3.09 
days, Msini = 0.93Mj, ap = 0.045 AU). The enhancement is clearly planet-related because 
it is not in phase with the stellar rotation rate Prot = 8 — 9 days (Tinney et al. 2001). Because 
the Ca II HK enhancement is seen only on the hemisphere facing the planet, a magnetic 
interaction is preferred over a tidal one. Since a unipolar inductor model (viz., the lo- Jupiter 
interaction; Zarka et al. 2001) does not fit the data as well (Saar et al. 2004), this appears to 
be the first observational evidence for an exoplanetary magnetosphere. These observations 
also provide an opportunity to probe the high-energy particle environment near EGPs since 
the interaction strength depends on the magnitudes of the stellar and planetary magnetic 



^ The current state of this rapidly advancing field is summarized by the International Astronomical 



Union's Working Group on Extrasolar Planets at http : //www . ciw . edu/boss/IAU/div3/wgesp/| and at the 



Exoplanet Encyclopedia at http://exoplajiet.eu/ 
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fields 5* and Bp (Saar et al. 2004). Note however that these emissions appear to be not 
phase locked, as is expected if the stellar magnetic fields are varying (Mclvor, Jardine, & 
Holzwarth 2006; Cranmer & Saar 2006). 

The peak strength of the HK emission fiux enhancement found is ~4% (Shkolnik et al. 
2004). Such low- level enhancements are difficult to detect and study in any detail. However, 
if the emission scales like other manifestations of stellar activity, coronal enhancements should 
be much larger: e.g., Ayres et al. (1995) found surface fiuxes Fx/F^oi oc (-Pciv/-^6o0^- Follow- 
up observations of HD 179949 in X-rays (Chandra/ ACTS; Saar et al. 2006) show significant 
spectral and temporal variability phased to the planetary orbit, but there is some residual 
ambiguity due to the poorly constrained stellar rotational period (Prot ~ 7 — 10 d) - some 
fraction of the variation may be due to changes in the underlying stellar activity. 

This picture has been reinforced by recent observations of similar activity on v And and 
a confirmation of the activity synchronization on HD 179949 (Shkolnik et al. 2008). There 
are also indications of the infiuence of a planet on stellar activity in r Boo and HD 189733, 
though the results are inconclusive. The newer observations also show that the possible 
infiuence of the giant planet on stellar activity is complex, intermittent, and prone to phase 
shifts (Saar et al. 2006, Shkolnik et al. 2008, Lanza 2008). 

In ^ we detail the X-ray data available for stars with giant plants. In ^ we carry 
out statistical searches for trends in the data as a function of ap, and show evidence of a 
significant deviation between extremal subsamples. In these analyses, we take into account 
the large numbers of censored data (that is, undetected stars) and properly include their 
effect. We discuss these results in §l]with careful attention to the biases in the sample. We 
summarize our results in ^ 



2. Data 

To datep over 230 star systems have been identified as possessing planets; various meth- 
ods such as spectroscopically detecting the wobble due to planetary orbits in radial velocity 
measurements (e.g., Butler et al. 1996) and detecting photometric dips in the light curve due 
to disk transit (e.g., Konacki et al. 2003,2004) have been used in these identifications. For 
the sake of simplicity and homogeneity (see §4.11) . we shall use for our sample only those stars 
which have had planets discovered or verified using the spectroscopic method. These stars 
are listed in Table [TJ In this list, we do not include those stars in which low-mass stellar or 
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brown-dwarf companions have been detected using the above methodso We have searched for 
X-ray counterparts of these stars in archival data from ASCA, EXOSAT, Einstein, ROSAT, 
XMM-Newton, and Chandra missions. Using the HEASARC BrowscI, we find matches in 
the ASCA Medium Sensitivity Survey (ascagis; Ueda et aL 2001), the EXOSAT resuhs 
(sc_cma_view; White & Peacock 1988), the Einstein 2-a catalog (twosigma; Moran et al. 
1996), the Einstein/IFC source catalog (einstein2e; Harris et al. 1994), the ROSAT/ERl 
complete results archive (roshritotal; Voges et al. 2001), the Brera Multi-scale Wavelet 
ROSAT/RRl catalog (bmwhricat; Panzera et al. 2003), the ROSAT/FSFC complete results 
archive (rospspctotal; Voges et al. 2001), the WGACAT (wgacat; White, Giommi, & An- 
gelini 1994), the ROSAT All-Sky Survey (RASS) Bright Sources catalog (rassbsc; Voges 
et al. 1996), the RASS Faint Sources catalog (rassf sc; Voges et al. 1999), the RASS A- 
K Dwarfs/Subgiants catalog (rassdwarf ; Hiinsch, Schmitt, & Voges 1998), and the XMM 
Serendipitous Source catalog (xmmssc; Pye et al. 2006). No matches have been found in the 
ChaMP database (Kim et al. 2004), and with the exception of e Eri, none of these stars have 
been observed with Chandra. 



^ While such a study is intrinsically interesting on its own merits, it is beyond the scope of our analysis. 
Binary star systems tend to be more X-ray active than single stars (e.g., Pye et al. 1994), possibly due to 
the higher angular momentum available or due to the manner of the evolution of the system, whereas the 
primary mechanisms for activity enhancements for stars with close-in EGPs appears to be tidal or magnetic 
disruption. Here we concentrate on a well-defined and limited sample of ostensibly single stars to explore 
the possible changes in their properties. 



* http : //heasarc . gsf c . nasa . gov/db-perl/W3Browse/w3browse . pi 
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Table 1. Stars with Giant Planets 



Name 


Spectral 


rny 


d 


ap- 


X-ray Flux*^ 


logio Lx 


eomments'^ 




Type 




[pc] 


[AU] 


[IQ-'^'-^ ergs cm-2 g-l 


[ergs s-i] 




HD 41004 B 


M2 


6.67 


43.0 


0.018 


11.6 ± 1.55 


29.41 ±0.06 


[138] P 


OGLE-TR-56 


G 


6.20 


1500 


0.023 


< 1.79 


< 30.26 


[57] P 


TrES-3 


G 


5.52 


292. 


0.023 


< 0.646 


< 29.82 


[91] P 


OGLE-TR-113 


K 


20.15 


1500 


0.023 


< 1.80 


< 30.45 


[58] P 


WASP-4 


G7V 


4.60 


300. 


0.023 


< 0.646 


< 29.84 


[135] P 


WASP-5 


G4V 


5.10 


297. 


0.027 


< 0.646 


< 29.83 


[1]P 


GJ 436 


M2.5 


5.49 


10.2 


0.029 


1.17 ± 0.525 


27.16 ±0.19 


[17] P 


SWEEPS- 11 


F9 


5.95 


2000 


0.030 


64.9 ± 4.55 


33.49 ±0.03 


[102] P 


OGLE-TR-132 


F 


5.00 


1500 


0.031 


< 1.78 


< 30.45 


[7]P 


WASP-2 


KIV 


9.16 


207. 


0.031 


< 0.646 


< 29.52 


[29] P 


HD 189733 


K1-K2 


10.08 


19.3 


0.031 


6.08 ± 0.183 


28.43 ±0.01 


[9] P E H X 


WASP-3 


F7V 


20.10 


223. 


0.032 


< 0.646 


< 29.58 


[96] P 


HD 212301 


F8 V 


13.60 


52.7 


0.036 


< 2.13 


< 28.85 


[63] P 


HD 63454 


K4 V 


12.57 


35.8 


0.036 


0.655 ± 0.285 


28.00 ±0.19 


[81] P 


TrES-2 


GOV 


3.73 


220. 


0.037 


< 1.46 


< 29.93 


[90] P 


XO-2 


KOV 


6.70 


149. 


0.037 


< 0.646 


< 29.23 


[10] P 


HD 73256 


G8/K0 


13.00 


36.5 


0.037 


2.15 ± 0.905 


28.53 ±0.18 


[129] P 


HAT-P-7 


F3 


10.68 


320. 


0.038 


< 0.646 


< 29.90 


[93] P 


WASP-1 


F7V 


9.36 


499. 


0.038 


< 0.646 


< 30.28 


[29] P 


HAT-P-3 


K 


11.40 


140. 


0.039 


< 0.646 


< 29.18 


[124] P 


HD 86081 


F8V 


3.23 


91.0 


0.039 


< 1.46 


< 29.16 


[50] P 


GJ 674 


M2.5 


10.42 


4.54 


0.039 


13.9 ± 2.05 


27.53 ±0.06 


[68] P N 


TrES-1 


KOV 


10.56 


157. 


0.039 


< 1.14 


< 29.53 


[105] P 


HD 83443 


KOV 


6.17 


43.5 


0.041 


< 1.70 


< 28.59 


[76] P 


HAT-P-5 


K 


10.17 


340. 


0.041 


< 0.646 


< 29.95 


[3]P 


Gl 581 


M3 


10.40 


6.26 


0.041/0.073/0.25 


< 1.64 


< 26.89 


[6]P 


HD 46375 


Kl IV 


8.71 


33.4 


0.041 


0.212 ± 0.036 


27.45 ± 0.07 


[72] N 


TWHya 


K8V 


11.86 


54.0 


0.041 


42.7 ± 0.147 


30.17 ±0.00 


[115] N P A 


OGLE-TR-10 


Gor K 


11.20 


1500 


0.042 


0.081 ± 0.047 


30.34 ±0.25 


[8]N 


HD 187123 


G5 


12.00 


50.0 


0.042 


< 0.073 


< 27.34 


[13] N 


HD 330075 


G5 


10.50 


50.2 


0.043 


0.047 ± 0.032 


27.15 ±0.30 


[75] N 


HD 149026 


GO IV 


10.50 


78.9 


0.043 


< 1.15 


< 28.93 


[108] P 


HD 2638 


G5 


8.42 


53.7 


0.044 


< 0.987 


< 28.53 


[81] P 


HAT-P-4 


F 


8.21 


310. 


0.045 


< 0.646 


< 29.87 


[60] P 


HD 209458 


GO V 


7.47 


47.0 


0.045 


0.039 ± 0.018 


27.02 ±0.20 


[22] N 


HD 179949 


F8 V 


8.06 


27.0 


0.045 


4.68 ± 1.12 


28.61 ±0.10 


[120] P 


T Boo 


F7 V 


5.80 


15.0 


0.046 


25.1 ± 0.0886 


28.83 ±0.00 


[12] N P E H A 


HD 75289 


GOV 


7.36 


28.9 


0.046 


< 1.68 


< 28.23 


[125] P 


BD-10 3166 


G4 V 


5.52 


100. 


0.046 


3.14 ± 1.01 


29.58 ±0.14 


[14] P 


Lupus-TR-3 


KIV 


6.29 


1780 


0.046 


< 0.646 


< 30.07 


[134] P 


OGLE-TR-111 


G8 


8.02 


1500 


0.047 


< 0.646 


< 29.69 


[105] P 


HD 88133 


G5 IV 


7.00 


74.5 


0.047 


< 1.08 


< 28.86 


[40] P 


XO-3 


F5V 


8.74 


260. 


0.048 


< 0.646 


< 29.72 


[49] P 


XO-1 


GIV 


8.20 


200. 


0.049 


< 0.646 


< 29.49 


[78] P 


TrES-4 


F 


7.61 


440. 


0.049 


< 1.50 


< 30.54 


[67] P 



-6- 



Table 1 — Continued 



Name 


Spectral 


my 


d 


ap'' 


X-ray Flux*^ 


logio Lx 


comments'^ 




Type 




[pc] 


[AU] 


[10-1'^ ergs cm-2 g-l 


[ergs s-l] 




HD 76700 


G6 V 


8.80 


59.7 


0.049 


< 2.20 


< 28.97 


[122] P 


HD 102195 


KOV 


6.69 


29.0 


0.049 


1.06 ± 0.423 


28.03 ±0.17 


[45] P 


OGLE-TR-182 


G2 


7.30 


2618 


0.051 


< 1.76 


< 30.63 


[97] P 


OGLE-TR-211 


F4 


7.57 


2177 


0.051 


< 1.83 


< 31.28 


[117] P 


HD 219828 


GOIV 


7.54 


81.1 


0.052 


< 0.989 


< 28.89 


[79] P 


51 Peg 


G2 IV 


7.30 


14.7 


0.052 


0.247 ± 0.052 


26.80 ±0.09 


[74] P 


HAT-P-6 


F 


7.18 


200. 


0.052 


< 0.646 


< 29.49 


[89] P 


HD 149143 


GOIV 


8.07 


63.0 


0.053 


< 1.47 


< 28.84 


[39] P 


SWEEPS-04 


F5 


6.42 


2000 


0.055 


< 0.646 


< 29.61 


[102] P 


HAT-P-1 


GOV 


4.70 


139. 


0.055 


< 0.646 


< 29.17 


[68] P 


HD 49674 


G5 V 


7.60 


40.7 


0.058 


< 1.59 


< 28.50 


[16] P 


V And 


F8 V 


5.48 


13.5 


0.059/0.83/2.5 


5.98 ± 0.920 


28.11 ±0.07 


[12] P 


HD 109749 


G3 IV 


8.33 


59.0 


0.064 


< 1.73 


< 28.86 


[39] P 


HD 168746 


G5 


7.44 


43.1 


0.065 


< 1.55 


< 28.54 


[94] P 


HAT-P-2 


F8 


7.51 


135. 


0.068 


< 0.646 


< 29.15 


[68] P 


HIP 14810 


G5 


8.05 


52.9 


0.069/0.41 


< 1.82 


< 28.78 


[18] P 


HD 118203 


KO 


7.57 


88.6 


0.070 


< 1.81 


< 29.23 


[139] P 


HD 68988 


GO 


8.45 


58.0 


0.071 


< 1.77 


< 28.85 


[132] P 


HD 162020 


K2 V 


6.45 


31.3 


0.072 


12.3 ± 2.58 


29.16 ±0.09 


[127] P 


HD 285968 


M2.5V 


6.94 


9.40 


0.073 


2.87 ± 0.889 


27.48 ±0.13 


[33] P 


HD 217107 


G8IV 


7.25 


37.0 


0.073/4.4 


0.203 ± 0.075 


27.52 ±0.16 


[34] N 


HD 185269 


GOIV 


5.70 


47.0 


0.077 


< 2.17 


< 28.76 


[51] P 


HD 69830 


KOV 


7.69 


12.6 


0.079/0.19/0.63 


1.56 ± 0.535 


27.47 ±0.15 


[66] P 


HD 130322 


KOV 


7.34 


30.0 


0.088 


0.361 ± 0.050 


27.59 ±0.06 


[125] N 


HD 108147 


F8/G0 V 


5.38 


38.6 


0.104 


0.356 ± 0.167 


27.80 ±0.20 


[94] N 


G186 


KIV 


8.16 


11.0 


0.110 


< 0.646 


< 26.97 


[98] P 


HD 4308 


G5 V 


7.90 


21.9 


0.114 


< 2.06 


< 28.07 


[130] P 


55 Cnc 


G8 V 


7.03 


13.4 


0.115/0.24/5.8/0.038 


0.493 ± 0.126 


27.02 ±0.11 


[12] E 


HD 27894 


K2 V 


6.74 


42.4 


0.122 


<4.51 


< 28.99 


[81] P 


HD 99492 


K2V 


7.25 


18.0 


0.123 


0.937 ± 0.374 


27.56 ±0.17 


[71] P 


HD 38529 


G4IV 


7.27 


42.4 


0.129/3.7 


4.96 ± 0.709 


29.03 ±0.06 


[35] P 


HD 195019 


G3 IV-V 


6.51 


20.0 


0.139 


0.052 ± 0.020 


26.40 ±0.17 


[34] N 


HD 192263 


K2 V 


7.24 


19.9 


0.150 


2.26 ± 1.06 


28.03 ±0.20 


[103] P 


HD 6434 


G3IV 


5.15 


40.3 


0.150 


< 3.10 


< 28.78 


[76] P 


HD 102117 


G6V 


6.78 


42.0 


0.153 


< 1.98 


< 28.62 


[65] P 


HD 17156 


GO 


7.28 


78.2 


0.159 


< 2.06 


< 29.18 


[46] P 


HD 33283 


G3V 


9.18 


86.0 


0.168 


< 1.99 


< 29.24 


[50] P 


Gliesc 876 


M4 V 


6.99 


4.72 


0.208/0.13/0.021 


0.837 ± 0.165 


26.35 ± 0.09 


[30] P X 


pCrB 


GOV or G2V 


5.95 


17.4 


0.220 


< 1.36 


< 27.69 


[12] P 


HD 11964 


G5 


6.92 


34.0 


0.229/3.2 


< 1.53 


< 28.32 


[18] P 


HD 224693 


G2IV 


7.95 


94.0 


0.233 


< 2.41 


< 29.41 


[50] P 


HD 43691 


GOIVV = 8.03 


5.91 


93.2 


0.240 


< 1.58 


< 29.21 


[46] P 


HD 37605 


KOV 


8.22 


42.9 


0.250 


< 2.67 


< 28.77 


[27] P 


HD 107148 


G5 


7.74 


51.3 


0.269 


< 0.547 


< 28.24 


[18] P 


HD 117618 


G2V 


8.31 


38.0 


0.280 


< 1.33 


< 28.36 


[123] P 
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Table 1 — Continued 



Name 



Spectral 
Type 



mv 



d 

[pc] 



[AU] 



X-ray Flux*^ 
[10~^^ ergs cm~^ s" 



logio Lx 

[ergs s^-"-] 



comments 



HD 3651 


KO V 


8.17 


11.0 


0.284 


1.23 ± 0.232 


27.25 ±0.08 


[38] H 


HD 74156 


GO 


8.03 


64.6 


0.294/3.8/1.0 


< 1.30 


< 28.81 


[84] P 


HD 219449 


Koni 


7.18 


45.0 


0.300 


< 0.619 


< 28.18 


[80] P 


HD 114762 


F9V 


7.98 


39.5 


0.300 


0.376 ± 0.049 


27.85 ±0.06 


[62] N 


HD 168443 


G5 


6.25 


37.9 


0.300/2.9 


< 1.56 


< 28.43 


[127] P 


HD 101930 


Kl V 


7.86 


30.5 


0.302 


< 1.94 


< 28.33 


[65] P 


HD 121504 


G2 V 


6.68 


44.4 


0.320 


< 1.41 


< 28.52 


[76] P 


HD 178911 B 


G5 


7.22 


46.7 


0.320 


1.66 ± 0.633 


28.64 ±0.17 


[137] P 


HD 16141 


G5IV 


7.83 


35.9 


0.350 


<2.46 


< 28.58 


[72] P 


HD 80606 


G5 


8.22 


58.4 


0.439 


< 1.80 


< 28.86 


[82] P 


HD 216770 


Kl V 


7.68 


38.0 


0.460 


< 1.96 


< 28.53 


[76] P 


HD 93083 


K3 V 


7.31 


28.9 


0.477 


< 1.34 


< 28.13 


[65] P 


70 Vir 


G4 V 


5.71 


22.0 


0.480 


0.458 ± 0.134 


27.42 ±0.13 


[69] H 


GJ 3021 


G6 V 


7.70 


17.6 


0.490 


28.3 ± 3.63 


29.02 ±0.06 


[83] P N 


HD 52265 


GOV 


7.79 


28.0 


0.490 


<2.05 


< 28.28 


[14] P 


HD 208487 


G2V 


6.45 


45.0 


0.490 


<4.35 


< 29.02 


[123] P 


HD 37124 


G4 V 


6.91 


33.0 


0.530/3.2/1.6 


< 1.70 


< 28.34 


[131] P 


HD 231701 


F8V 


7.60 


108. 


0.556 


<2.10 


< 29.47 


[41] P 


HD 155358 


GO 


6.40 


42.7 


0.628/1.2 


< 1.69 


< 28.57 


[28] P 


HD 73526 


G6 V 


5.06 


99.0 


0.660/1.0 


< 1.56 


< 29.26 


[122] P 


ksi Aql 


G9nib 


8.08 


62.7 


0.680 


< 123. 


< 30.76 


[110] P 


HD 75898 


GO 


6.41 


80.6 


0.737 


< 1.35 


< 29.02 


[101] P 


HD 8574 


F8 


9.01 


44.2 


0.760 


<3.75 


< 28.94 


[95] P 


HD 134987 


G5 V 


7.38 


25.0 


0.780 


< 1.31 


< 27.99 


[131] P 


HD 104985 


G9 III 


7.48 


102. 


0.780 


< 1.47 


< 29.26 


[107] P 


HD 81688 


Koni-iv 


7.65 


88.3 


0.810 


<2.74 


< 29.41 


[110] P 


HD 169830 


F8 V 


6.63 


36.3 


0.810/3.6 


0.887 ± 0.053 


28.15 ±0.03 


[76] P 


HD 40979 


F8 V 


5.94 


33.3 


0.811 


< 1.64 


< 28.34 


[36] P 


HD 150706 


GO 


7.77 


27.2 


0.820 


8.57 ± 0.714 


28.88 ±0.04 


[126] P I 


HD 202206 


G6 V 


6.80 


46.3 


0.830/2.5 


< 1.61 


< 28.62 


[127] P 


HD 12661 


G6 V 


6.03 


37.2 


0.830/2.6 


< 1.10 


< 28.26 


[35] P 


4 Uma 


Kim 


6.06 


62.4 


0.870 


1.47 ± 0.101 


28.84 ±0.03 


[31] H P 


HD 89744 


F7 V 


8.10 


40.0 


0.890 


0.597 ± 0.079 


28.06 ±0.06 


[59] P 


HR810 


GOV pecul. 


6.18 


15.5 


0.910 


19.1 ± 2.73 


28.74 ±0.06 


[61] P 


HD 59686 


K2ni 


4.21 


92.0 


0.911 


< 1.86 


< 29.28 


[80] P 


GJ 317 


M3.5 


8.04 


9.17 


0.950 


<3.20 


< 27.25 


[52] P 


HD 92788 


G5 


7.83 


32.8 


0.970 


< 0.963 


< 28.09 


[35] P 


HD 142 


Gl IV 


7.70 


20.6 


0.980 


<3.13 


< 28.20 


[121] P 


HD 156846 


GOV 


8.23 


49.0 


0.990 


< 1.91 


< 28.74 


[119] P 


HD 177830 


KO 


7.10 


59.0 


1.000 


< 1.77 


< 28.87 


[131] P 


ChaHaS 


M6.5 


8.58 


160. 


1.000 


53.7 ± 2.02 


31.22 ±0.02 


[48] A N 


HD 122430 


K3III 


8.99 


135. 


1.020 


< 1.61 


< 29.54 


[111] P 


HD 28185 


G5 


7.24 


39.4 


1.030 


< 2.21 


< 28.61 


[104] P 


HD 175541 


G8IV 


9.80 


128. 


1.030 


<2.95 


< 29.76 


[46] P 


HD 100777 


KO 


4.44 


52.8 


1.030 


< 0.504 


< 28.23 


[86] P 



Table 1 — Continued 



Name 


Speetral 


my 


d 


ap^ 


X-ray Flux'' 


logio Lx 


comments'^ 




Type 




[pc] 


[AU] 


[10-i« ergs cm-2 s-^ 


[ergs s-i] 




HD 142415 


Gl V 


9.42 


34.2 


1.050 


2.88 ± 0.577 


28.61 ±0.09 


[76] P 


HD 108874 


G5 


7.81 


68.5 


1.051/2.7 


< 0.603 


< 28.53 


[16] P 


HD 4203 


G5 


9.97 


77.5 


1.090 


0.629 ± 0.294 


28.66 ±0.20 


[132] E 


HD 128311 


KO 


8.41 


16.6 


1.099/1.8 


10.2 ± 1.35 


28.52 ±0.06 


[16] P H U 


HD 33564 


F6 V 


9.00 


21.0 


1.100 


1.30 ± 0.455 


27.84 ±0.15 


[43] P 


HD 210277 


GO 


8.05 


21.3 


1.100 


< 1.32 


< 27.85 


[70] P 


HD 99109 


KO 


5.10 


60.5 


1.105 


< 0.548 


< 28.38 


[18] P 


HD 192699 


G8IV 


5.80 


67.0 


1.160 


<4.36 


< 29.37 


[68] P 


HD 210702 


KlHI 


7.68 


56.0 


1.170 


< 1.45 


< 28.73 


[68] P 


HD 27442 


K2 IV a 


8.69 


18.1 


1.180 


0.826 ± 0.304 


27.51 ±0.16 


[15] P 


HD 82943 


GO 


5.94 


27.5 


1.190/0.75 


< 1.13 


< 28.01 


[76] P 


HD 188015 


G5IV 


5.67 


52.6 


1.190 


< 2.46 


< 28.91 


[71] P 


HD 125612 


G3V 


6.73 


52.8 


1.200 


< 1.48 


< 28.69 


[41] P 


HD 114783 


KO 


8.65 


22.0 


1.200 


< 0.861 


< 27.70 


[132] P 


HD 154857 


G5V 


8.65 


68.5 


1.200 


< 2.20 


< 29.09 


[77] P 


HD 221287 


F7 V 


7.91 


52.9 


1.250 


3.35 ± 1.25 


29.05 ±0.16 


[68] P 


HD 20367 


GO 


8.69 


27.0 


1.250 


20.8 ± 0.165 


29.26 ±0.00 


[126] N P 


HD 147513 


G3/G5V 


7.79 


12.9 


1.260 


51.5 ± 1.16 


29.01 ±0.01 


[76] P X 


HIP 75458 


K2III 


6.54 


31.5 


1.275 


0.228 ± 0.065 


27.43 ±0.12 


[42] P 


HD 4113 


G5V 


8.03 


44.0 


1.280 


< 5.14 


< 29.08 


[119] P 


HD 171028 


GO 


7.88 


90.0 


1.290 


<2.63 


< 29.41 


[106] P 


HD 17092 


Koin 


7.84 


109. 


1.290 


< 1.64 


< 29.37 


[88] P 


HD 19994 


F8 V 


5.26 


22.4 


1.300 


<2.25 


< 28.13 


[76] P 


HD 167042 


KlHI 


8.10 


50.0 


1.300 


< 1.73 


< 28.71 


[46] P 


HD 41004 A 


Kl V 


7.22 


42.5 


1.310 


11.6 ± 1.55 


29.40 ± 0.06 


[128] P 


HD 222582 


G5 


6.86 


42.0 


1.350 


< 0.773 


< 28.21 


[131] P 


HD 20782 


G2 V 


6.30 


36.0 


1.360 


< 1.68 


< 28.42 


[56] P 


HD 65216 


G5 V 


8.05 


34.3 


1.370 


< 0.646 


< 27.96 


[76] P 


HD 160691 


G3 IV-V 


5.45 


15.3 


1.500/4.2/0.090/0.92 


0.423 ± 0.061 


27.07 ±0.06 


[15] P 


HD 141937 


G2/G3 V 


1.15 


33.5 


1.520 


< 1.38 


< 28.27 


[127] P 


HD 183263 


G2IV 


9.40 


53.0 


1.520 


<3.71 


< 29.10 


[71] P 


HD 47536 


Kl HI 


7.72 


121. 


1.610/5.2 


< 1.83 


< 29.51 


[112] P 


HD 114386 


K3 V 


7.98 


28.0 


1.620 


< 1.12 


< 28.02 


[76] P 


HD 23079 


F8/G0 V 


8.00 


34.8 


1.650 


<3.86 


< 28.75 


[121] P 


HD 4208 


G5 V 


8.21 


33.9 


1.670 


< 1.51 


< 28.32 


[132] P 


16 Cyg B 


G2.5 V 


5.95 


21.4 


1.680 


< 1.46 


< 27.90 


[25] P 


HD 62509 


KOIIIb 


8.70 


10.3 


1.690 


0.687 ± 0.022 


26.94 ±0.01 


[100] N P E A 


V391 Peg 


sdB 


7.18 


1400 


1.700 


< 1.55 


< 31.56 


[lie'' ] P 


HD 5319 


G5IV 


7.48 


100. 


1.750 


< 1.60 


< 29.28 


[101] P 


HD 70573 


G1-1.5V 


8.08 


45.7 


1.760 


2.57 ± 0.423 


28.81 ±0.07 


[114] E P 


OGLE-05-071L 


M5 


9.00 


2900 


1.800 


1.05 ± 0.189 


32.03 ±0.08 


[68^ ] P 


HD 13189 


K2II 


7.61 


185. 


1.850 


< 1.59 


< 29.81 


[47] P 


HD 45350 


G5IV 


6.36 


49.0 


1.920 


< 1.60 


< 28.66 


[71] P 


eps Tau 


KO III 


8.04 


45.0 


1.930 


0.615 ± 0.175 


28.17 ±0.12 


[109] P 


HD 11977 


G8.5 HI 


8.13 


66.5 


1.930 


<2.83 


< 29.17 


[113] P 
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Table 1 — Continued 



Name 


Spectral 


mv 


d 


ap'' 


X-ray Flux*^ 


logio Lx 


comments'^ 




Type 




[pc] 


[AU] 


[10-1^ ergs cm-2 g-i 


[ergs s-i] 




HD 81040 


G2/G3 


8.93 


32.6 


1.940 


1.21 ± 0.508 


28.19 ±0.18 


[118] P 


HD 111232 


G8V 


7.72 


29.0 


1.970 


< 2.18 


< 28.34 


[76] P 


HD 132406 


GOV 


5.41 


71.0 


1.980 


< 1.91 


< 29.06 


[46] P 


HD 159868 


G5V 


6.54 


52.7 


2.000 


<2.11 


< 28.84 


[92] P 


HD 213240 


G4 IV 


8.24 


40.8 


2.030 


< 1.84 


< 28.56 


[104] P 


Gamma Cephei 


K2 V 


7.11 


11.8 


2.044 


0.550 ± 0.068 


26.96 ±0.05 


[26] P 


HD 187085 


GOV 


8.74 


45.0 


2.050 


< 2.14 


< 28.71 


[56] P 


HD 16175 


GO 


8.06 


59.8 


2.070 


< 1.58 


< 28.83 


[46] P 


HD 190647 


G5 


7.01 


54.2 


2.070 


< 2.01 


< 28.85 


[68] P 


HD 114729 


G3 V 


5.74 


35.0 


2.080 


< 1.24 


< 28.26 


[16] P 


NGC 2423 3 


KHI 


7.10 


766. 


2.100 


< 2.89 


< 31.31 


[64] P 


OGLE-05-390L 


M4 


8.33 


6500 


2.100 


< 0.646 


< 30.82 


[68^ ] P 


HD 10647 


F8V 


8.80 


17.3 


2.100 


4.51 ± 0.260 


28.21 ±0.02 


[128] P 


HD 164922 


KOV 


7.57 


21.9 


2.110 


< 1.54 


< 27.95 


[18] P 


47 Uma 


GOV 


8.51 


14.0 


2.110/7.7 


0.583 ± 0.050 


27.13 ±0.04 


[11] N 


HD 10697 


G5 IV 


3.31 


32.6 


2.130 


1.30 ± 0.520 


28.22 ±0.17 


[131] P 


HD 2039 


G2/G3 IV-V 


5.40 


89.8 


2.190 


<2.68 


< 29.41 


[122] P 


HD 170469 


G5IV 


17.40 


65.0 


2.240 


< 677. 


< 30.87 


[41] P 


OGLE-06-109L 


K9 


10.04 


1490 


2.300/4.6 


< 0.646 


< 28.81 


[44d ] p 


HD 136118 


F9 V 


7.40 


52.3 


2.300 


< 1.71 


< 28.75 


[37] P 


HD 190228 


G5IV 


23.00 


66.1 


2.310 


< 2.52 


< 29.12 


[95] P 


HD 11506 


GOV 


22.30 


53.8 


2.350 


< 1.70 


< 28.77 


[41] P 


Gj 849 


M3.5 


18.34 


8.80 


2.350 


1.91 ± 0.486 


27.25 ±0.11 


[19] E 


HD 50554 


F8 


20.17 


31.0 


2.380 


< 0.076 


< 26.94 


[95] N 


NGC 4349 No 127 


Kill 


15.78 


2176 


2.380 


< 0.646 


< 31.56 


[64] P 


HD 23127 


G2V 


17.98 


89.1 


2.400 


< 4.88 


< 29.67 


[92] P 


HD 196050 


G3 V 


16.08 


46.9 


2.500 


< 2.25 


< 28.77 


[76] P 


18 Del 


G6III 


16.08 


73.1 


2.600 


4.18 ± 0.836 


29.43 ± 0.09 


[110] P 


HD 106252 


GO 


16.84 


37.4 


2.610 


< 0.287 


< 27.68 


[95] P 


HD 196885 


F8V 


14.80 


33.0 


2.630 


< 1.83 


< 28.38 


[18] P 


HD 216435 


GOV 


16.56 


33.3 


2.700 


1.08 ± 0.290 


28.16 ±0.12 


[55] N 


kappa CrB 


KllVa 


21.20 


31.1 


2.700 


< 1.29 


< 28.17 


[53] P 


HD 216437 


G4 IV-V 


99.00 


26.5 


2.700 


0.179 ± 0.055 


27.18 ±0.13 


[76] N 


HD 23596 


F8 


21.30 


52.0 


2.720 


< 1.59 


< 28.71 


[95] P 


14 Her 


KOV 


17.40 


18.1 


2.770 


< 0.528 


< 27.32 


[73] N 


OGLE-05-169L 


MO 


18.80 


2700 


2.800 


< 0.646 


< 28.47 


[68'* ] P 


HD 142022 A 


KOV 


19.83 


35.9 


2.800 


< 1.93 


< 28.47 


[32] P 


HD 66428 


G5 


11.10 


55.0 


3.180 


< 1.90 


< 28.84 


[18] P 


HD 39091 


Gl IV 


4.50 


20.5 


3.290 


0.423 ± 0.050 


27.33 ± 0.05 


[54] P 


HD 70642 


G5 IV-V 


4.50 


29.0 


3.300 


0.058 ± 0.019 


26.77 ±0.14 


[21] N 


e Eri 


K2 V 


4.50 


3.20 


3.390 


108. ± 0.360 


28.12 ±0.00 


[20] N P E H X A U 


HD 117207 


G8VI/V 


12.40 


33.0 


3.780 


< 1.21 


< 28.20 


[71] P 


HD 30177 


G8 V 


4.50 


55.0 


3.860 


< 3.20 


< 29.06 


[122] P 


HD 50499 


G IV 


4.50 


47.3 


3.860 


< 3.70 


< 28.99 


[133] P 


HD 190360 


G6 IV 


4.09 


15.9 


3.920/0.13 


< 0.286 


< 26.94 


[85] N 
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We also carried out a targeted survey of some of the stars with planets using XMM- 
Newton (PI: V.Kashyap). The stars were chosen from the extreme ends of the distribution of 
ap, to provide a contrast between stars with close-in planets and stars with distant planets 
(see §3.31) ■ These stars are listed in Table [21 We used the XMM Science Analysis System 
(SAS vT.O.O; 20060628_1801) to reduce the data and obtained source counts within circular 
cells of radius 20" and a background estimated from nearby regions devoid of sources. In the 
cases where no excess X-ray emission above the background was detected at a significance 
of > 0.997 (corresponding to a Gaussian-equivalent 3a detection), count rate upper limits 
were calculated as described by Pease et al. (2006). 

Overall, we find that 70 of the stars in the sample have been detected serendipitously 
or in pointed observations (Table [1]). For those stars left undetected, we first determine 
an upper limit from the RASS data, by estimating the number of counts required for a 3cr 
detection (see Pease et al. 2006) given the observed count rate in the 0.1-2.4 keV band at 
the location of the source and the accumulated exposure time (in the all-sky maps from the 
survey; Snowden et al. 1997). If a pointed XMM-Newton observation exists, we use an upper 
limit derived from that observation. We further impose a limit of -r^^ < 10~^ since stars do 
not exceed this limit on average (though such limits may be exceeded on occasion when a 
large flare occurs). For the detected stars, we compute a nominal counts-to-energy conversion 
factor, cecf assuming a coronal source spectrum with similar temperature components of 2 
and 5 MK and an absorption column of 10^^ cm~^. For XMM-Newton data, we computed the 
cecf at a higher temperature (10 MK) to ensure that similar values of fluxes are obtained for 
sources detected with both the MOS and the pn detectors. Adopting a single value of cecf 
for all coronal sources observed with a given detector introduces systematic errors of ~ 30%, 
comparable to the statistical errors present in the measurements, and significantly smaller 
than the intrinsic variability in the sources and other systematic errors present in the the X- 
ray luminosity functions; the dominant source of uncertainty in our results is sample variance 
arising from the intrinsic variations in the X-ray luminosity functions. Using WebPIMMS 
v3.4 and XSPEC vlO, we estimate that the cecf for 

ASCA/GIS = 3 X 10"^°, 

Einstein/IPC = 1.8 x 10"", 

EXOSAT/LE = 6.5 x 10"^°, 

ROSAT/BRl = 2.8 x 10"", and 

ROSAT/FSFC = 6.5 x 10"^^ ergs cm"^ ct"^ 

XMM/EPIC - pn = 1.9 x 10"^^ ergs cm'^ ct"^ 



|http : //heasarc . gsf c . nasa.gov/Tools/w3piimns . html 
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Table 1 — Continued 



Name 



Spectral mv 
Type 



d 

[pc] 



[AU] 



[10- 



X-ray Flux'' 
^'^ ergs crn^'^ s 



logio Lx 

[ergs s^-'-] 



comments 



HD 89307 


GOV 


4.09 


33.0 


4.150 


< 0.883 


< 28.06 


[18] P 


HD 72659 


GOV 


11.79 


51.4 


4.160 


< 1.92 


< 28.78 


[16] P 


HD 154345 


G8V 


11.98 


18.1 


4.190 


< 1.49 


< 27.76 


[68] P 


SCR 1845 


M8.5 V 


10.64 


3.85 


4.500 


3.20 ± 0.254 


26.75 ±0.03 


[5d]P 


OGLE235-MOA53 


K5 


12.60 


6025 


5.100 


< 0.646 


< 29.61 


[68^ ] P 


2M1207 


M8 


12.26 


52.4 


46.000 


< 0.646 


< 25.91 


[23'' ] P 


GQ Lup 


K7eV 


11.30 


140. 


103.00 


2.53 ± 0.397 


29.77 ±0.07 


[87d ] P 


ABPic 


K2 V 


11.18 


45.6 


275.00 


39.8 ± 1.69 


30.00 ± 0.02 


[24d ] P 


UScoCTIO 108 


M7 


9.80 


145. 


670.00 


< 1.40 


< 29.55 


[2^1 P 



N E U 



^Semi-major axes of planetary orbits 

''The X-ray flux adopted from the best available measurement (see text) 

'^The reference for the planetary detection (sec below), and the X-ray mission in which the star was detected; A=ASCA/GIS, 
X=EXOSAT/LE, H=ROSAT/HRI, E=Einstcin/IPC, P=ROSAT/PSPC, N=XMM-JVewton 

''Planets around stars have not been detected via the Radial Velocity method 
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et al. (2007); [117] Southworth (2008); [118] Sozzetti et al. (2006); [119] Tamuz et al. (2007); [120] Tinney et al. (2000); [121] Tinney et 
al. (2001); [122] Tinney et al. (2002); [123] Tinney et al. (2005); [124] Torres et al. (2007); [125] Udry et al. (2000); [126] Udry et al. 
(2002a); [127] Udry et al. (2002b); [128] Udry et al. (2003a); [129] Udry et al. (2003b); [130] Udry et al. (2005); [131] Vogt et al. (2000); 
[132] Vogt et al. (2002); [133] Vogt et al. (2005); [134] Weldrake et al. (2008); [135] Wilson et al. (2008); [136] Wolszczan & Frail (1992); 
[137] Zucker et al. (2002); [138] Zucker et al. (2004); [139] da Silva et al. (2005). 
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XMM/EFIC - MOS = 8.4 x 10"'^ ergs cm-^ ct 



in the 0.1-4.5 keV passbandS in case a star has been detected multiple times with different 
instruments, the flux measurement with the best S/N is reported in column 8 of Table [H 
Similarly, if observed multiple times but never detected, we report the lowest of the upper 
limits computed for that star. 



^ There are slight differences in the conversion factors among different catalogs from the same instrument, 
e.g., rassbsc and wgacat. However, these discrepancies are much smaller than the uncertainty caused by 
the source spectral distributions, and has no effect on our analyses. 
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Table 2. XMM-Newton Data 



Name 


ap 


ObsID/Rcvolution 




Exposure [s] 




Count rate [ct ks" 


-'] 




[AU] 




MOS-1 


MOS-2 


PN 


MOS-1 


MOS-2 


PN 


HD 46375 


0.041 


0304202501/1071 


8215.9 


8464.2 


6123.0 


2.13 ±0.77 


< 1.5 


11.2 ±1.9 


HD 187123 


0.042 


0304203301/1166 


14586.7 


15368.0 


9382.7 


< 1.25 


<1.11 


<3.25 


HD 330075 


0.043 


0304200401/1037 


15309.5 


15772.6 


13400.2 


<0.8 


< 0.8 


< 1.7 


HD 217107 


0.073 


0304200801/0995 


7536.3 


7797.8 


5497.1 


< 1.1 


<1.14 


4.1 ±1.5 


HD 130322 


0.088 


0304200901/1028 


6295.6 


6758.4 


4208.7 


4.4 ±1.1 


4.14 ±1.1 


19 ±2.7 


HD 190360 


0.128 


0304201101/0985 


3078.0 


3837.3 


2573.9 


<9.76 


<2.93 


< 13.9 


HD 195019 


0.138 


0304201001/1167 


10327.4 


10259.4 


8211.5 


< 1 


< 1 


2.73 ±1.2 


47UMa 


2.11 


0304203401/1191 


8564.7 


7735.4 


6131.2 


7.81 ±1.1 


8.1 ±1.2 


30.7 ±2.6 


HD 50554 


2.38 


0304203201/1163 


1125.7 


1117.0 


1098.1 


<6.3 


< 4.4 


< 7.9 


HD 216435 


2.70 


0304201501/1166 


2045.0 


2793.5 


183.05 


12.8 ±3.5 


8.21 ±2.4 


84.4 ±30.4 


HD 216437 


2.70 


0304201601/0979 


4163.5 


4396.5 


2677.1 


< 2.3 


<2.15 


9.4 ±2.9 


14 Her 


2.77 


0304202301/1059 


4565.7 


4570.9 


- 


< 5.66 


<7.83 


- 


HD 70642 


3.30 


0304201301/1159 


13509.7 


13489.1 


10763.8 


<0.9 


< 0.9 


3.1 ±0.99 


HD 33636 


3.56 


0304202701/1054 


10807.2 


10847.9 


6825.0 


4.1 ±0.84 


4.2 ±0.85 


19.4 ±2.1 
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3. Analysis 

Here we carry out a series of statistical tests on the sample of stars with close-in EGPs 
that were first detected, or at least verified, using the radial- velocity method. We first 
analyze the sample as a whole ( §3.ip . then search for correlations with orbital radius ap 
( §3.2p . and finally search for differences within extremal subsamples to increase the contrast 
of the signal ( §3.3p . 



3.1. Sample Statistics 

Of the >230 stars in our sample thus far identified as possessing EGPs, 70 are found to 
be X-ray emitters. If we exclude giants, we are left with a sample of 180 main sequence star 
systems with 46 X-ray detections. This is a smaller fraction of X-ray bright stars than is 
found for field stars in the solar neighborhood (e.g., Schmitt 1997 finds 95% of F stars and 
83% of G stars within 13 pc are detected with ROSAT, while Maggio et al. 1997 find that 
~50% of G stars within 25 pc that were observed with Einstein are detected). 

A direct comparison of these surveys with our star sample is however misleading; one 
must account for the intrinsic variations in X-ray luminosities and the inhomogeneous scatter 
in distances in our sample. To investigate this, we carry out a Monte Carlo simulation on our 
sample of stars, fixing each star at its given distance, but allowing their X-ray luminosities 
to vary. These luminosities are obtained by adopting the known X-ray luminosity functions 
(derived from statistically complete samples of nearby field stars using Einstein data; they 
thus include all the variations known to exist for coronal sources; cf. Kashyap et al. 1992). 
We generate 1000 realizations of X-ray fiuxes at Earth using these luminosity functions 
for each star, and consider them to be detected in X-rays if the fiux at Earth is > 2 x 
]^Q-i3 gj^gg Qj2i~^ s"^ (corresponding to the typical RASS sensitivity; Schmitt 1997). The 
resulting frequency distributions of imputed detections can then be compared to the observed 
frequency distribution of X-ray detections in our sample. This comparison is shown in 
Figure [H where the frequency distributions of actual X-ray detections in our sample (dash- 
dotted line) and of simulated X-ray detections for a nominal set of field stars located at the 
same distances as the stars in our sample (dashed line, with error bars derived from averaging 
over the Monte Carlo realizations) are shown as a function of the planetary orbit size. Also 
shown are the fractional residuals between the actual and predicted number counts (solid 
histogram), which shows a weak trend towards a higher efficiency of detection for stars with 
close-in planets (dashed line; slope=— 0.1 ± 0.03). A quantitative comparison between the 
two distributions yields a reduced x^ of ~ l-S; denoting that they are marginally statistically 
distinguishable. 
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Fig. 1. — Distribution of observed and expected number of stars, as a function of the semi- 
major axis, ap of the closest planet. Left; The solid histogram represents the sample of main 
sequence stars among those listed in Table [H while the dot-dashed histogram represents the 
subset of those stars which have been detected in X-rays. The dashed histogram with 
error-bars denoted by vertical line segments is the distribution generated by adopting known 
X-ray luminosity functions and carrying out Monte Carlo simulations for stars at the specific 
distances given in Table [H Right: The histogram shows the fractional difference between 
the predicted and the actual number of X-ray detections in the sample. The solid horizontal 
line depicts the case of a perfect match, while the dashed line shows the best-fit linear fit to 
the residuals. This shows that there may exist a weak trend in the data towards a higher 
number of X-ray detections for stars that have close-in planets. 
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A comparison of the X-ray luminosity distribution of the sample of stars with planets, 
compared to a nominally unbiased distribution from field stars is shown in Figure [2l The 
expected X-ray distributions are shown as the shaded region, and the actual luminosity 
distribution (derived as a Kaplan- Meier estimator; see e.g., Schmitt 1985, Feigelson & Nelson 
1985) as the solid histogram. A formal Kolmogorov-Smirnoff test between the two shows that 
the null hypothesis that the data are drawn from the full sample cannot be ruled out. Most 
of the differences can be attributed to detections over a range of Lx ~ 3 x 102"^"^^ ergs s~^ 
While not definitive, this again suggests the need for a more sensitive analysis. 



3.2. Correlation with ap 

We show above ( §3.ip that, taken as a whole, the set of main sequence stars with EGPs is 
similar in gross characteristics to, though not identical to, the field star sample. Given that, 
we next consider possible variations within our sample, to test whether the EGPs have any 
measurable effect on levels of X-ray emission. Because frequency distributions are heavily 
model dependent and are limited to integer values, they do not have sufficient statistical 
power to detect weak variations in the properties of the stars. We expect that the effect of 
EGPs will increase as the orbital distance decreases, and vice versa; we therefore check for 
direct correlations of various sample parameters with the orbital distance of the EGP. 

We also note that the sample is not volume limited, and that this may introduce a 
number of complications into the analysis. We therefore consider an additional filter to 
obtain a subset of the main-sequence stars that are within 60 pc. The number distribution 
of these stars is uniform within this distance (see Figure [3]), and this subsample comes closest 
to a statistically complete volume limited sample. 

We expect a priori that tidal and magnetic effects due to close-in giant planets will 
manifest themselves as a trend in stellar activity as a function of the planetary orbit ap. 
We have thus carried out detailed statistical analyses to measure correlations of the X-ray 
luminosity Lx and the activity indicator j^ with ap. We have carried out these tests 
for the full sample of main-sequence star systems which have had planets detected via the 
Radial Velocity method, as well as for the subsample of stars which are known to be X-ray 
emitters, and for the subsample of stars which lie within 60 pc. If a given star has multiple 
planets detected, we choose the planet with the smallest semi-major axis as defining ap. 
We compute both Pearson's and Spearman's Rank Correlation coefficients, and report the 
results in Table [31 The errors on Pearson's coefficient are derived via bootstrapping by 
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.0 27.5 28.0 28.5 29.0 29.5 50.0 

logio(Lx [ergs s"^]) 



Fig. 2. — X-ray luminosity distribution function. The fraction of stars with a given luminos- 
ity Lx is shown for both the main sequence sample from Table [1] (solid histogram) as well 
as for the expected distribution based on a model of the Galaxy (shaded region). The width 
of the shading represents the statistical error present, both in terms of the numbers of stars 
as well as the uncertainties inherent in the construction of the X-ray luminosity functions. 
The two cannot be statistically distinguished under a Kolmogorov-Smirnoff test. 
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Fig. 3. — The distribution of the number of main-sequence star systems with distance. The 
cumulative number within a given distance is plotted as the sohd curve. The vertical dashed 
line represents the distance limit of the statistically complete subsample at 60 pc, within 
which the star systems are uniformly distributed. 
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sampling with replacement |j The significance of Spearman's p is denoted by the p-value, 
which measures the probabihty that the observed value of p can be obtained as a chance 
fluctuation. 

We show the distribution of X-ray detections and upper limits for the full sample in 
Figure m as a function of ap. No trend is visually discernible here in Lx{ap)- Detailed 
correlation analysis (see Table [3]) shows that there may be a slight negative correlation: for 
the full sample, Spearman's Rank correlation coefficient is p = —0.09 ± 0.05 (the la error 
bar on p is computed via 10000 Monte Carlo simulations with the uncertainties in Lx taken 
into account at each iteration). This is weakly significant, but the correlation becomes less 
significant when smaller samples are considered (for the X-ray detections, p = 0.01, with 
p = 0.96 even when uncertainties in Lx are ignored, and for the volume limited sample, 
p = —0.08 with p = 0.38). We also compute Pearson's correlation coefficient, which uses the 
information of the actual values of the luminosities and not just their rank order, and thus 
produces statistically more powerful estimates. As before, we obtain error bars by carrying 
out Monte Carlo simulations that take into account the uncertainties in Lx', we adopt a 
Gaussian sampling function with a equal to the measured error for detected sources and 
a fiat distribution in logLx for undetected sources. We find that Pearson's coefficient is 
-0.01 ± 0.06 for the full sample, -0.03 ±0.1 for the X-ray detected sample, and 0.02 ± 0.3 
for the volume limited sample. 

A similar weak correlation is found when the correlation of the activity indicator -r^^ 
with ap is considered. For the set of all stars, Pearson's coefficient is — 0.04±0.11 (Spearman's 
p = —0.25 ± 0.001), for only X-ray detected stars, it is 0.16 ±0.1 (Spearman's p = 0.01; 
p = 0.96), and for the volume limited sample, 0.26 ±0.12 (Spearman's p = —0.08; p = 0.38). 
We consider the activity indicator in greater detail below (see §4.11) in establishing the 
magnitude of the sample bias. 

Since we do not expect that Lx and ap are linearly related, and because the dynamic 
range is large, it is not surprising that the correlation tests are contradictory and inconclusive. 
A more sensitive analysis is therefore required, and we adopt a technique wherein the contrast 
in the data is maximized (see §3.31) by increasing the lever arm between extremal subsamples. 

In addition to the X-ray luminosity Lx, we also consider the correlations of the stellar 
distance and stellar radius with ap. We expect a negative correlation to exist with distance 
to the star, simply because it is easier to detect planets and activity around nearby stars. 



^ Where measurement errors are available (e.g., for Lx), wc include them in the simulations by sampling 
with a Gaussian of standard deviation equal to the measurement error and a mean equal to the maximum 
likelihood value; upper limits are dealt with using a uniform distribution. 
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Fig. 4. — X-ray luminosities of the sample stars as a function of the orbital distance of the 
closest EGP. The downward arrows mark upper limits based on non-detections in RASS. 
The X-ray detections are marked with an 'x'; the errors on the measured fluxes are denoted 
by vertical bars. 
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Table 3. Correlation of orbital semi-major axis ap with stellar parameters'* 

Parameter Dataset Pearson's Coefficient*^ Spearman's p'^ 

Lx Full sample -0.01 ± 0.06 -0.09 ± 0.05 

X-ray detections -0.03 ± 0.1 0.01 {p = 0.96) 

volume limited 0.02 ± 0.29 -0.08 {p = 0.38) 

logio(;^) Full sample -0.04 ±0.1 -0.09 ± 0.05 

X-ray detections 0.16 ± 0.09 0.006 {p = 0.96) 

volume limited 0.26 ± 0.12 -0.08 {p = 0.38) 

d^ Full sample -0.19 ±0.03 -0.2 (p = 0.003) 

X-ray detections -0.28 ± 0.1 -0.07 [p = 0.6) 

volume limited 0.09 ± 0.06 -0.07 {p = 0.5) 

R^ Full sample -0.023 ± 0.045 -0.05 {p = 0.47) 

X-ray detections -0.05 ± 0.09 0.10 {p = 0.48) 

volume limited -0.05 ± 0.1 0.17 {p = 0.06) 



^Stars from programs such as OGLE, SWEEPS, TrES, WASP, etc., which 
initially detected planets via the photometric transit method and not the 
radial velocity method, have been excluded from these tests. 

^Pearson's linear correlation coefficient; errors derived via bootstrapping 
by sampling with replacement. Where possible, measurement errors and 
upper limits are accounted for via Monte Carlo simulations. 

'^Rank correlation of two populations; the p value denotes the two-sided 
significance of its deviation from by random chance, i.e., small values 
indicate significant correlation. Where p-values are not quoted, the la 
error on the correlation coefficient, computed via Monte Carlo simulations, 
is shown. 
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This is borne out by the correlation analysis: we find that Pearson's coefficient is — 0.19±0.03 
(Spearman's p = —0.21, with p = 0.003) for the entire sample, and for the X-ray detected 
sample alone, Pearson's coefficient is —0.28 ±0.1 (Spearman's p = —0.07, p = 0.6). For 
the volume limited sample, in contrast, we obtain a Pearson's coefficient of 0.09 ± 0.06 
(Spearman's p = 0.07, p = 0.46). Note that here we determine the errors via bootstrapping 
by sampling with replacement. We expect to find no correlation with stellar radius, and 
this expectation is borne out by the correlation analysis. For the full sample, Pearson's 
coefficient is —0.02 ±0.04 (Spearman's p = —0.05, p = 0.47); for the X-ray detected sample, 
-0.05 ± 0.09 (p = 0.1, p = 0.48); and for the volume limited sample, -0.05 ± 0.1 (p = 0.17, 
p = 0.06). 



3.3. Extremal Subsamples 

In order to clarify the variations in the X-ray activity statistically, we consider the 
extreme cases of stars with planets that are close to the primary (the "close-in" sample) and 
compare their luminosities with those of stars with planets at large distances (the "distant" 
sample) o We choose these subsamples from the volume limited sample, i.e., main sequence 
star systems which are within 60 pc. This signal is obviously dependent on the contrast 
between the close-in and the distant samples, and in order to obtain the best contrast it is 
necessary to choose samples that are as far apart in their range of ap as possible and yet 
contain as large a sample of stars as possible. A comprehensive investigation of the best 
such separation is not feasible for this sample because of the large number of X-ray non- 
detections and the variety of tests that we carry out on the data. However, for our purposes it 
is sufficient to determine whether there exists subsamples which show the requisite contrast, 
at some separations, even if it is not necessarily the optimal one. (Note that we do consider 
below the effect of varying the sizes of the samples.) We therefore adopt an ad hoc separation 
based on the population distribution (Figured]): we choose ap < 0.15 AU (corresponding 
to a dip in the frequency distribution of stars as a function of orbital distance ap) for the 
close-in sample, which results in a sample size of 40 stars, 20 of which are detected in X-rays. 
By limiting the lower bound of the distant sample such that there is separation of an order 
of magnitude difference in ap between the limits of the two samples, we set for the distant 
sample, ap > 1.5 AU; 8 of 38 stars in this subsample are detected in X-rays. This range 
ap has the additional advantage that similar numbers of stars are found in each set. The 



^We are precluded from directly comparing our sample with that of an independent sample of field stars 
because at this stage it is not possible to be certain that the putative field star sample is devoid of close-in 
EGPs. 
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mean X-ray luminosities for these two subsamples are found to be well separated, with the 
close-in sample being significantly X-ray brighter (see Figure [5]). We comment on alternative 
choices further below (see §4.21) and explicitly show (Figure EJ that this specific choice does 
not affect our results. 

Because of the large number of upper limits in the dataset, we cannot carry out simple 
hypothesis tests to verify whether these two samples are derived from the same parent 
distribution. Instead, we carry out Monte Carlo realizations of the sample as before ( §3.2p . 
using a Gaussian error distribution with measured errors for the detected stars and a flat 
uniform distribution in the log-scale for the undetected stars. We compute sample means for 
each set of realizations, and the distribution of these means allows us to determine whether 
the two .ample, are similar or different! 

We find that the mean X-ray luminositjo of the close-in sample is < Lx{cip < 0.15) >= 
;Lq28.49±o.o9 g^gg g-1^ ^^^ ^Yi^^ ^f |-]^g distant sample is < Lx{ap > 1.5) >= lO^^-ssio.is g^gg g-i_ 

The mean level of emission for the close-in sample is significantly greater than for the distant 
sample, and as shown in Figure El the two samples are found to be different at > 99% 
confidence level o 

As discussed above, our choice of the limiting ap for the two subsamples is ad-hoc. We 
have therefore considered the effect on our result of of varying the subsample bounds of ap. 
The mean Lx for various subsamples is obtained using Monte Carlo simulations as above, 
for various ranges of ap, and are shown in Figure O Two types of subsamples are considered: 
one that includes all stars with planetary orbital radii < ap^i^g^_in, from which the mean X- 
ray luminosity < Lx{ap < apdose-in) > is obtained (upper shaded curve), and another that 
includes all stars with orbital radii > ap^isianti which results in < Lx{ap > ap^istant) > (lower 
shaded curve). As can be seen, the sample that includes close-in planets is invariably more 
X-ray intense than the sample that includes distant planets for almost all possible choices 



^ Note that because fewer than half the stars in the sample are detected in X-rays, we are precluded from 
using sample medians as a summary estimator of the subsamples; the median is not a robust estimator in 
this case. 

^°Here and henceforth the enclosing brackets "<>" denote the mean value of a quantity; note that within 
these brackets we often have numerical conditional expressions such as "<" (less than) and ">" (greater 
than), enclosed within parentheses. 

^^ Note that the nearby strong X-ray source eEri is in the distant subsample. A planet with a period 
~ 10 yr was detected around it by Campbell, Walker, & Yang (1998), but this detection remains contro- 
versial due to the large intrinsic scatter of ~ 20 m s^^ in the radial velocity curves (see Marcy et al. 2002; 



http://exoplcinets.org/esp/epseri/epseri.shtinll. Excluding eEri would decrease the mean Lx of the 



distant sample even further, increasing the separation between the two distributions. 
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Fig. 5. — Mean Lx for the extreme samples. The distribution of simulated means for the 
set of stars with the closest giant planet having semi-major axis ap < 0.15 AU is shown as 
the solid histogram, and its counterpart for stars that have giant planets with ap > 1.5 AU 
is shown as the dashed histogram. The overlap between the two distributions is < 1%. 
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of < ap close-in and > apiiistant- Note that these results do not change in any quahtative way 
if the dM stars (which contribute significantly to the high mean Lx point at small ap) are 
excluded from the sample. 

We caution that this observed difference in the mean Lx between stars with close-in 
and distant planets cannot be entirely attributed to the effect of close-in EGPs; there are 
selection biases inherent in the sample that must be accounted for. By considering indirect 
sample ensemble properties, we argue below (see §4.ip that the bias inherent in our selected 
sample has a small effect. Based on physical grounds (see § ^4.3p . we expect that close-in 
giant planets could have a significant effect on the X-ray activity level of the primary, and 
the trend seen in the activity level as a function of the size of the planetary orbit (semi-major 
axis ap) must in large part be due to the effect of the close-in giant planet. 



4. Discussion 

4.1. Sample Bias 

The sample of candidate stars for which planet searches are conducted is subject to some 
subtle biases. Some of these biases have the effect of masking the signature of planet induced 
activity and are difficult to quantify since the planet detection processes are numerous, the 
programs are still incomplete, and the rates of false positives and false negatives are unknown. 
We may however determine the approximate extent of these biases by studying the ensemble 
properties of the sample of stars with detected planets. Here we describe these biases, and 
consider their effect on our ability to detect intrinsic trends in X-ray activity. Furthermore, 
because the set of EGPs is dominated by those identified with the radial velocity method, 
we shall limit our sample to those stars which have been verified to have planetary systems 
by this method, and will thus consider only those biases introduced by that method. 

1. Spectral Homogeneity Most of the stars for which planets are searched for are solar 
like. This is advantageous to our study since the stars considered here are relatively 
homogeneous; therefore effects that arise due to changes in the convective turnover 
timescales leading to changes in the nature of the magnetic dynamo at the high- and 
low-mass ends of the coronal sequence, or due to the changing evolutionary states of 
the systems may generally be ignored. We further homogenize the sample by limiting it 
to main sequence stars within 60 pc (additional types of filtering to further homogenize 
the samples has no effect on our results; see §4.21) . 

Because the sample of stars is homogeneous, we do not expect any correlations between 
the stellar radii i?* and ap to be present, and indeed we find that the data are consistent 
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Fig. 6. — Comparison of mean Lx for extremal subsamples. The width of the shading 
represents the la statistical error on the mean Lx derived via bootstrapping and Monte 
Carlo simulations (see text). The upper shaded region represents subsamples that include 
close-in planets, < Lx{ap < o-p dose-in) > and the lower shaded region, the subsamples that 
include distant planets < Lx{ap > ap distant) >■ The mean values in both cases are shown 
as thick dark lines marked with diamonds within the shaded regions. The distributions in 
Figure [5] correspond to vertical cuts across the shaded regions here, at apdose-in = 0.15 (upper 
region, for < Lx{ap < 0.15) >) and apj^^gig^^t = 1-5 (lower region, for < Lx{ap > 1-5) >). 
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with these expectations (Table [3]). We have also carried out a full analysis of the 
extremal subsamples for various subsets of the full dataset (see Figure [10]) and show 
that our results are robust to such selections. We thus conclude that in our sample 
there is no bias present due to stellar type or size effects. 

2. Stellar Distance For a given planetary mass, the detectability of planets decreases 
as the distance of the planet from the star {ap) increases, since the radial velocity 
amplitude decreases. Since detectability in general decreases with distance to the star 
{d^,), we expect to find more stars with small ap at larger d^,. In the sample of stars 
with detected planets, we therefore expect that ap and d* are anti-correlated. Indeed, 
we find a slight negative correlation between ap and d^, though it is weak (see Table [3]). 
This bias exists simply as a result of the limitations of measurement statistics and is 
independent of any stellar activity effects. Thus, if we assume that X-ray properties of 
stars are independent of the distance limit, the only effect on the X-ray data will be to 
have a larger fraction of higher upper limits at smaller ap (due to the larger numbers 
of more distant stars) o Thus, stellar distance should have no effect on the ensemble 
properties of the X-ray luminosities of stars from such a sample. In any case, we carry 
out most of our analysis with the volume limited sample, which is explicitly designed 
to remove any such effect. 

3. Selecting for reduced activity The process of planet detection via the radial velocity 
method (RV; e.g., Butler, et al. 1996) is limited by the amount of intrinsic RV jitter 
that may be caused by stellar activity (Saar & Donahue 1997, Saar et al. 1998; see also 
Baluev 2008), and hence planets are preferentially detected when they are close-in, 
massive, and the primary star is not active. 

Because radial velocity signatures are hard to detect around active stars, there is a 
tendency to select candidate stars for planet search against activity. Thus, a priori 
there is an expectation that the sample of stars with detected planets would have 
lower activity levels than field stars of the same type. But as shown in §3.H the X-ray 
emission from these stars is by and large consistent with field star emission levels. In 
any case, in order to avoid introducing accidental biases by comparing the sample of 
stars with EGPs to those without detected EGPs (note that a lack of detection does 
not preclude the existence of planets around the star - even if the star is part of the 
planet search program, there may be a planet that is less massive or more distant than 



^^ In our correlation analyses we have tested for the effect that these higher upper limits may have by 
using uniform distributions (bounded above by the value of the upper limit) in both log space and normal 
space to describe the censored data and find no qualitative difference between the two cases. 
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the sensitivity that has been reached), we have chosen to compare the extrema of a 
single distribution, viz., the sample of stars known to possess giant planets. Thus any 
global selection effect that may exist in the sample towards lesser activity will apply 
equally to both subsamples and its effects are irrelevant for this study. 

Also note that in some types of planet detection methods, persistent stellar activity 
fixed to an active longitude can mimic the signature of a planet (Lanza et al. 2008). 
This results in false detections of planets at the stellar orbital period. However, this 
has an effect on photocentric methods and have no effect on radial velocity or transit 
methods currently in use. 

4. Intra-sample trend in inhibition of activity If stellar activity inhibits planet 
detection, then it will inhibit it selectively. Planets at large ap are preferentially 
detected around stars with weak activity, since high activity may mask the RV signal 
of the distant planet. Thus, in a sample of stars selected based on the existence of 
planets, those stars with distant planets are presumptively less active. This bias thus 
produces the same signature in activity trends that we search for (see §3.3p . and thus 
interferes with our study. 

How can this bias be quantified? Consider that the cause of the bias is the excess jitter 
in velocity induced by magnetic activity, which is related to the energy deposited in the 
corona, which in turn is tracked by the ratico of X-ray to bolometric luminosity -j^^^. 
This jitter serves to mask the RV signal, which is due to wobble from gravitational 
effects and is therefore inversely proportional to the orbital semi-major axis of the 
planet. Thus, the bias will manifest itself as a positive correlation between the formally 
independent parameters -j^^^ and -^. Note that this correlation is not necessarily linear 
because it is dependent on the observation process, the cumulative observation time, 
etc. Any existing correlation also does not imply causality, and indeed the velocities 
that are generated by these two processes are quite different. Nevertheless, we expect 
that stars with large values of -jj^- will primarily be present in our sample provided they 
are also accompanied by small values of ap. Thus, a measurement of the correlation 
between these two variables serves as a measurement of the sample bias. For the sample 
of main sequence stars within 60 pc that have EGPs detected via the RV method, we 
find that 

logio T^ = (-5.65 ± 0.04) + (0.22 ± 0.05) ■ log^o — , (2) 

^Bol CLp 

where ap is in units of [AU]. Here the errors are derived via Monte Carlo bootstrapping 
on the X-ray flux measurement errors. Thus, we find that there does exist a weak. 



13 



■j^ is preferred for this bias measurement because it is surface area independent. 
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but statistically significant, correlation between j^ and -^. This trend is shown in 
Figure [7l Conservatively assuming that all of this correlation is attributable to the 
inherent bias, we can then estimate its effect on the observed activity trend. That is, 
for each of the stars in the close-in and distant samples (in §3.3p . given ap and LboI 
we can compute a predicted Lx{ap,LBoi) that can then be used in the place of the 
measured Lx to carry out the same analysis. The result of this calculation is shown 
in Figure [HI which shows that the close-in and distant samples differ by ^ la, and 
account for only a factor of ~ 2 of the observed difference between the two samples. 
It is important to note that the magnitude of this bias will vary for every subsample, 
and that the magnitude and direction of the bias will be case specific. 

Note that j^ is the preferred proxy for activity because it is surface area independent. 
However, because there are no R^: biases in the sample, Lx and j^ are similarly 
correlated with ap. Thus, a priori, we expect that accounting for the selective inhibition 
bias of Equation [2] in the ap vs Lx dataset should result in a complete removal of 
any luminosity difference between the extremal samples, i.e., that ^<^ d°tant>^ ^ ^' 
That this is manifestly not so (Figure [9]) is an indication that the differences in the 
distribution of Lx for the extremal samples is systematic, and an inherent property of 
the observed sample. 

Furthermore, the bias measurement is conservative, i.e., we include some of the desired 
signal within the estimate of the bias, and thus we overestimate the magnitude of the 
bias. Because this bias is removed from the observed ratios of the mean luminosities, 
our results provide a conservative estimate of the true magnitude of the effect of planet- 
induced activity. 

5. Intrinsic sample variance Other parameters known to affect the level of stellar 
activity are age, rotation, the Rossby number, metallicity, temporal variability, etc. 

Stellar ages are not known with sufficient accuracy for us to consider them here, and 
it is possible that there exists a trend where ap is correlated with stellar age. But 
note that such hidden variables do not affect the results derived here, and simply point 
to a more complex explanation of the connection between close-in planets and stellar 

activity. 

It is also as yet unclear what type of effect a massive Jupiter-type planet will have on 
these parameters. For instance. Rotational synchronization by tidal interactions may 
be thought to increase v sin i and thereby coronal activity for close-in EGPs compared 
to distant EGPs. However, not only is this as likely to spin-down a fast rotating star as 
to spin-up a slow-rotator, but also the timescales for the synchronization of the stellar 
rotation period with the planetary orbital period are very large to begin with, and 
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Fig. 7. — The distribution of log -^^^^ with log^. The X-ray detections are marked with 
'x's, with error bars denoted with vertical lines; the X-ray non-detections are shown with 
downward arrows. Also shown are the best-fit regression line (solid curve) and the envelope 
bound of the bootstrapped regression curves (dotted lines). Since more active stars have 
higher radial velocity jitter, we expect fewer planets to be detected around stars at high j^^^ 
unless they also have high — ; this correlation tracks an important sample bias (see text). 
This trend is removed from the rest of the analysis. 
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increase as ap^ for distant planets (see e.g., Drake et al. 1998 and references therein), 
implying that most stars in the sample are not synchronized to the orbits of their 
close-in giant planets. We therefore conclude that both theoretical and observational 
prejudice points to evenly scattered values for these stellar parameters in our sample. 



4.2. Activity Enhancement 

We use the parameterization of the bias described above (see §4.ip to determine the 
residual signal, i.e., the magnitude of the excess X-ray activity in the close-in subsample 
compared to the distant subsample that can be attributed to the effect of EGPs. In order 
to simplify the calculation, we consider the ratio of the average luminosities of the close-in 
and distant subsamples calculated during the Monte Carlo simulation described above (see 
analysis in §3.31) . The data from Figure [5] are shown as the solid histogram in Figure [9], and 



<7 
>2.8) 



confirm that the close-in subsample is on average more X-ray luminous by a factor 4 
where the bounds on the number indicate the extent of the asymmetrical Icr-equivalent 
errors. After correcting for the sample bias (by dividing the original ratio by the ratio of 
the luminosities as predicted by Equation [2]), we find that the bias-corrected luminosity 
ratio decreases, but is > 1 at a significance of 0.95. This constitutes a definite detection of 
activity enhancement, by a factor 2.3^^3 (dashed histogram in Figure [HI). This remaining 
enhancement can be attributed to the effect of close-in giant planets on their parent stars. 
We therefore conclude that the close-in and distant samples indeed show differing X-ray 
activity levels and that a factor ~ 2 can be attributed to the effect of close-in EGPs. Note 
that the existence of this residual enhancement does not in itself allow us to conclude that 
there is a causal connection between the closeness of giant planets to their primary stars and 
the activity levels on those stars. However, in conjunction with the Call HK enhancements 
observed by Shkolnik et al. 2003), our results suggest that such a connection may be present. 

As a further check on the sensitivity of our analysis to the samples used, we compute the 
same ratio as in Figure for a number of different subsamples. The resulting distributions of 
the ratios of mean luminosity for close-in versus distant subsamples are shown in Figure [TOl 
in all cases the close-in subsamples have larger < Lx > than the distant subsamples. 

Note that the above result is dependent on the bounds chosen for ap for both subsamples. 
Clearly, the contrast between the two extremal sets would be heightened if the range of ap 
defining the two subsamples is shrunk further, and vice versa. This must also be balanced 
by the decreased robustness of the result due to the consequent dearth of X-ray detections 
in the subsamples. A rigorous estimate of the "best" split between the close-in and the 
distant samples is not feasible because of the large number of stars that remain undetected 
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in X-rays and the need to compute the bias factor (Equation [2]) separately for each case. 
Nevertheless, as we have shown above (Figure [6]) there exists a persistent difference in the 
mean luminosity between the close-in and distant subsamples, and it favors the suggestion 
that stars with close-in planets tend to be more active. The best way to improve this result 
is to obtain X-ray detections of more stars at the extreme ranges of ap. 



4.3. Dependence of the X-ray excess on ap 

For tidally interacting close binaries in a sample that included RS CVn stars, Schrijver 
& Zwaan (1991) found that the X-ray surface flux Fx oc a/o^, where a* is the distance to 
the cooler component. Comparing the average X-ray luminosities for close-in and distant 
subsamples that we derive above (see § §3.3114.21) . after sample biases have been corrected for, 
we find that 

<i^x_> oc < ap >-°■3^±°■l^ (3) 

which is a smaller dependence than that found by Schrijver & Zwaan, but the large error 
bar prevents a definitive conclusion. 

It is well known that binary stars are generally more X-ray active than single stars 

of the same type and rotation rate (see, e.g., Zaqarashvili, Javakishvili, & Belvedere 2002 

and references therein). Numerous mechanisms have been proposed to account for this 

so-called "overactivity," generally based on tidal and magnetic interactions. These studies 

were extended to the case of stars with EGPs and brown dwarfs by Cuntz et al. (2000). 

They suggested that the energy generation due to tidal interactions is proportional to the 

gravitational perturbation 

A^ ^ Mp 2Rl 

g, M, {ap - R,y ^ '^' 

where M* is the mass of the star and Mp that of the EGP, and the height of the tidal bulge, 

9* 

Mp , o 
oc -i-Rtap-^ 

oc ap~^ . (4b) 

Saar et al. (2004) estimated the energy released via reconnection during an interaction of 
the planetary magnetosphere with the stellar magnetic field, 

Fint OC B^BpVrelO.p'"' 

OC ap'"" , (5) 



-33- 



where 5* and Bp are the stellar and planetary magnetic fields, Vrei is the relative velocity 
between them that produces the shear in the magnetic fields and leads to the reconnection. 
Here, n = 3 very close to the star and n = 2 farther away (in the "Parker spiral"). Thus, if 
the magnitude of the enhancement and the stellar magnetic field can be measured for these 
systems, then the planetary magnetic environment can be investigated. 

Since our measurements are an ensemble average, i.e., they do not take into account 
variations in mass, magnetic fields, and relative velocities, we cannot directly verify the 
above models. Nevertheless, we can obtain a rough estimate for the variation of the "excess" 
luminosity with ap by assuming that 

5Lx oc ap^ (6) 

where 6Lx is the difference between the actual luminosity and a luminosity unaffected by a 
close-in EGP. For the former, we use the values simulated during the Monte Carlo analysis 
of §3.3^ and approximate the latter with the average Lx derived for the distant sample. 
Fitting straight lines to log 6Lx vs logap, and accounting for the bias in the same manner 
as described above in §4.2[ we find that x ~ — 1. The error on this value is however very 
large (~ 100%). Therefore, we conclude that while the data are qualitatively consistent with 
the scenarios proposed by Cuntz et al. (2000), a reliable test of the theory is not feasible 
without more observations. 



5. Summary 

We have searched for X-ray emission from a sample of 230 stars with known giant planets 
with a view to characterizing the effect of the EGPs on the parent star. We first find that 
the overall sample of stars with known EGPs is similar in gross X-ray properties to field 
star samples (see §3.11) . and thus provides a representative sample of X-ray stars. We carry 
out a careful search for statistical trends with various parameters (see §3.2p and find some 
evidence that the activity levels of stars with close-in giant planets is higher than for stars 
with planets located further out, though the correlations are contradictory and inconclusive. 
We then carry out a more powerful test by analyzing in detail two extremal subsamples (see 
where we compare the X-ray emission from nearby (distance<60 pc) main-sequence 



stars with close-in EGPs against similar stars with much more distant EGPs. For the sake 
of definiteness, we choose a close-in sample where the orbital semi-major axis ap < 0.15 
AU, and a distant sample where ap > 1.5 AU (reducing the separation between the samples 
increases the number of stars considered, but reduces the contrast between the samples). We 
verify that our adopte ranges of ap are not special by varying the limiting ranges of ap over 
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which the subsamples are defined, and find that invariably the close-in samples have X-ray 
luminosities higher than that of the distant sample. 

The above result must be understood in the context of selection biases in our sample 
of stars. In §4.11 we demonstrate that observational biases account for about half of the 
observed differences seen in the data. After these biases are accounted for, we find that the 
close-in sample is more active by a factor of 2^^ 3 on average. This result holds even when 
the full data are filtered with different selection criteria. 

The robustness of the result is limited by the large number of systems yet to be detected 
in X-rays. New observations by Chandra or XMM-Newton would consequently improve the 
statistics and would constrain the magnitude of this effect with better precision. We note 
that a simple model of the X-ray emission enhancement suggests an interaction strength 
proportional to the product of the stellar and planetary magnetic fields, B^ Bp at their point 
of interaction. This suggests that if B^, is known or can be estimated, Bp for exoplanets can 
potentially be studied. Since the point of interaction between Bp and S* depends in part 
on the strength of the stellar wind, close-in wind properties can also potentially be probed 
by such observations. 

This research has made use of data obtained through the High Energy Astrophysics 
Science Archive Research Center Online Service, provided by the NASA/Goddard Space 
Flight Center, and of the SIMBAD database, operated at CDS, Strasbourg, France. We 
thank Dr. Jean Schneider and the Exoplanets Encyclopedia, and the Geneva ExtraSolar 
Planets group for their invaluable online compilation of EGP resources. We also thank 
Pete Ratzlaff for assistance with computing counts-to-energy conversion factors. This work 
was supported by NASA grant NNG05GJ63G for XMM GO support, NASA-AISRP grant 
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Fig. 8. — Same as Figure [5l but using Lx{ap,LBoi) estimated from the regression analysis 
on j^ with ^ (Equation [21). The two distributions differ by ~ 2a. 
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Fig. 9. — Distributions of < Lx(close — in) > / < Lx(distant) >. The distribution of the 
ratios of simulated Lx for the close-in {ap < 0.15 AU) and distant {ap > 1.5 AU) samples 
is shown as the solid histogram. The same ratios, modified to account for the sample bias 
(see §4.11) are shown with dashed lines. The means and standard deviations of the two 
distributions are denoted with vertical and horizontal lines respectively (solid and dashed, 
for the two distributions). The data show that the two extremal samples differ by a factor 
~ 4 (between 2.8 to 7 times), and after sample bias is accounted for, by a factor ^ 2 (between 
1.3 and 4 times). 
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Fig. 10. — Effect of variations of sample composition on tlie distribution of < 
Lx(close — in) > / < Lx(distant) >. Tlie ratios of < Lx{ap < 0.15) > / < Lx{ap > 1-5) > 
are compared for different representative subsets of tlie full sample of stars with detected 
planets. In all plots, the baseline comparison from Figure [H for the main-sequence sample 
within 60 pc, is shown as the dashed histogram. The solid histograms are: (a) all X-ray 
detections, (b) the full sample of all stars with EGPs detected via the RV method, (c) the 
nearby main-sequence sample excluding dM stars, and (d) the full sample as in (b), exclud- 
ing giants. In all cases, the bias correction (see §4.ip is done separately for each sample. In 
all cases, the evidence is strong that stars with close-in giant planets are more active than 
those with distant planets. 



